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Abstract 

We have perfomed studies of the K\ — ► 7r + 7r~7 Direct Emission (DE) and 
Inner Bremsstrahlung (IB) vertices, based on data collected by KTeV during 
the 1996 Fermilab fixed target run. We find ai /a 2 = -0.737 ± 0.034 GeV 2 
for the DE form-factor parameter in the /5-propagator parametrization, and 
report on fits of the form-factor to linear and quadratic functions as well. We 
concurrently measure T(K% -> vr+ vr" 7 , E* > 20 MeV)/r(i^£ -» tt+tt') = 
(20.8 ± 0.3) x 10~ 3 , and a K° L -» vr+vr^ DE/(DE + IB) branching ratio of 
0.683 ±0.011. 

PACS numbers: ll.30.Er, 13.25.Es, 13.40.Hq, 14.40.Aq 
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The decay K\ — > 7r + 7r~7 (Figure |l|) is a potential new window into the phenomenon 
of CP-violation This decay arises primarily from the CP-violating Electric Dipole 

(El) "Inner Bremsstrahlung" (IB) and the CP-conserving Magnetic Dipole (Ml) "Direct 
Emission" (DE) processes. KTeV recently reported j|] the observation of a CP-violating 
angular asymmetry between the 7t + tt~ and e + e~ decay planes in the K\ — > 7r + 7r~7*, 7* — > 
e + e~ mode, arising from the interference of DE and IB photon polarization states [gj-|7|]. If 
an El DE term were also present, this would generate a direct CP-violating effect. 

Understanding the precise photon energy (E*) spectrum from the DE amplitude is crucial 
to interpretation of CP-violating effects and may shed light on particular chiral models |§. 



Previous experiments PJlQfl have observed evidence for an energy shift in the DE E* spec- 
trum, interpreted as evidence for an i?*-dependent form-factor modification to the pure-Mi 
DE amplitude. Here we report on the first measurement of this form-factor from direct fits 
to the data, using the rare decay mode K\ — > 

We consider two separate form-factor parametrizations. Historically, the p-propagator 
form m 

F = 1 — o 9~r — 7^ + 0-2 (1) 

(m 2 p - m 2 K ) + 2m K E* 

has been used, where ai/a 2 is the quantity of interest. More generally the form factor may 
be expressed clS db Taylor series in E* 

1 t_e; s_e; 2 

M K + M 2 K 



^=(1 + T7 2 + T7f + --- ( 2 ) 



where the interesting quantitites are the slope r and the quadratic parameter s. Recent 
theoretical work [|TT1 , |i"2| has attempted to calculate these parameters. 



The data presented here was collected by KTeV operating in the E832 configuration fig 
during the 1996 Fermilab fixed target run. A proton beam of intensity ~ 3 x 10 12 protons 
per 19 second spill incident at an angle of 4.8 mr on a BeO target was employed to produce 
two nearly parallel K\ beams. In E832, one of these beams was incident on an active Kg 
regenerator. Data collected in the regenerator beam was ignored for this analysis, except 



as a check of our understanding of the IB decay amplitude. The configuration of the KTeV 
E832 spectrometer consisted of a vacuum decay region, a magnetic spectrometer with four 
drift chambers, photon vetoes, a trigger scintillator bank, a Csl electromagnetic calorimeter, 
and a muon veto bank. 

Signal K\ —> 7r + 7r~7 and normalization K Q L — > 7r + 7r~ candidates were selected from the 
two-charged-track trigger. Offline, the sample was further purified by requiring the presence 
of two reconstructed tracks with a good vertex-^ 2 within the fiducial aperture of the detector, 
and no significant activity in the photon veto counters. The tracks' kinematics were required 
to be inconsistent with A — > p + ir~ decays. The energy deposited by the charged pions in 
the calorimeter was required to be less than 0.85 x the spectrometer momentum in order to 
eliminate backgrounds from K e3 decays. Candidate events were required to have a decay 
vertex between 110 and 156 meters downstream of the target and a total energy between 20 
and 160 GeV. 

K\ — > 7T + 7r _ 7 candidates were subject to the additional criterion of requiring that the 
quantity 

P 2 = (Mr ~ Ml» - Mlf - 4M 2 M 2 , - AM^U ^ 

4((P 2 u+iie) u 

be negative, i.e. by requiring that the tt° momentum be imaginary under a K\ — > 7r + 7r~7r° 
hypothesis. This cut virtually eliminates the — > 7r + 7r~7r background to the K\ — > 7r + 7r _ 7 
event sample. In addition, at least one non-track-associated cluster in the calorimeter was 
required to possess > 1.5 GeV of energy, and to be at least 3 cm removed from the calorimeter 
edges. This "photon cluster" had to be at least 20 cm from the nearest track projection 
to reject background from K\ — > 7t + tt~ events accompanied by pion hadronic showers in 
the calorimeter. The photon was required to have an energy of at least 20 MeV in the 
three-body center of momentum. 

Figures ^| (a) and (b) show the vacuum beam K\ s —> 7r + 7r~7 data after final cuts. A 
very clean signal of 8,669 K\ s — > 7r + 7r~7 events is achieved, with a background of about 
0.5%. Approximately 0.6% of the events in the peak in Figure § are residual K$ — > 7r + 7r~7 
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decays from Kg's generated at the target. A total of 4,482,706 K\ — > 7r + 7r~ events were 
accumulated for the normalization sample with a 0.1% background. 

To extract the Direct Emission form-factor, we consider the distribution in E* the 
photon energy in the center of momentum (Figure |3]). We wish to extract the relative 
contributions of the Inner Bremsstrahlung and DE terms, as well as the energy shift of the 
DE spectrum due to the presence of a form factor. We assume that El IB and form-factor- 
modified Ml DE are the only significant contributions to the decay. The possibility of an 
interference term arising from the presence of an El DE contribution is also briefly discussed 
below. 

We perform a MINUIT Jl4| ^-minimization fit to combine Monte Carlo DE and IB E* 
distributions, and extract the relative DE and IB contributions to the data as well as the 
DE form-factor. In the fit, the form factor parameters (a 1 /a 2 or r and s) and the ratio 

/ = I™. (4) 

J (T DE + T IB ) [) 

are allowed to float simultaneously. A x 2 is formed from comparison of the resultant DE+IB 
summed Monte Carlo histogram with the data, and minimized to obtain the best fit result. 

Table | summarizes the numerical fitting results for the three form-factor parametriza- 
tions. For the Taylor series parametrization, fits were performed with the quadratic param- 
eter both fixed at zero (floating a single parameter, r{) and allowed to vary (floating two 
parameters, r 2 and s^)- Data-Monte Carlo agreement in the E* distribution is shown in 
Figure |3|, for the p-propagator form. Figure [| illustrates the expected effect of the various 
parametrizations on the pure Ml direct emission spectrum. 

We see in the data presented in Table | and Figure |] clear evidence for a modification to 
the pure-Mi DE spectrum. All fits are good, though the p-propagator hypothesis stands 
out slightly: It gives the best \ 2 for a single-parameter fit, and the two-parameter fit 
results are in good agreement with the values r 2 = —2.70 and s 2 = 3.87 obtained by Taylor- 
expanding the p-propagator form. The size of this energy shift, in particular the need to 
take into account terms of second-order in E* is not currently understood within the chiral 



perturbation theory model for this decay [[U] . 



Additional contributions to the photon energy spectrum are expected from CP-violating 
higher-order multipole contributions to the Direct Emission amplitude. One possible con- 
sequence of these multipole terms is the presence of a charge-asymmetry in the ir + versus 
7r~ Dalitz plot We exclude asymmetries larger than 2.4% at 90% c.l. with the present 
data. 

One might also expect a contribution to the E* spectrum from CP-violating El DE, 
which could interfere with the El IB. We have allowed for such a term, by performing a 
separate fit in which we assumed an El DE amplitude constant in E* and searched for the 
corresponding interference contribution [|l!| to the E* spectrum. The form-factor parameter 



and / are allowed to float simultaneously, while fixing the IB rate at its theoretical value 
of 7.00 x 10~ 3 Based on this fit, we set an upper limit of T in /{Tde + ^ib + ^in) < 
0.30 (90% c.l.) on the contribution to the decay rate of the CP-violating interference term. 
Fin/Fall = 0.30 corresponds to a 22% decrease in the T DE /T\ LL ratio. 

Background subtraction under the E* distribution (with shape determined by study of 
sideband) resulted in no statistically significant change in either the form-factor parameter 
or the DE/(DE + IB) ratio. Accidental activity in the detector was also found to have 
no statistically significant effect. The form factor parameters are found to have a slight 
(1.9%) sensitivity to variations of the lower E* cutoff (nominally 20 MeV) used in the fit. 
We also assign systematic uncertainties of 1.8% due to detector acceptance, and 1.4% due 
to the effects of uncertainties in the calorimeter photon energy scale. Bias due to Level 3 
trigger inefficiency in the 1996 charged data set ||il||T5| is not significant in this analysis at 
the present level of statistics. 

Note (Table |) that the measured DE/(DE + IB) ratio is insensitive to the particular 
choice of form-factor parametrization. Uncertainty in the detector acceptance contributes a 
0.7% systematic uncertainty to this ratio, while variations in analysis cuts contribute 0.6%. 
The systematic uncertainty due to the calorimeter photon energy scale is 0.2%. 

We find the absolute K® — > 7r + 7r~7 branching ratio by normalizing the signal to the 
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K\ — > 7r + 7r _ channel. We determine the ratio of K\ — > 7T + 7r~/Ki — > 7r + 7r~7 acceptances 
using a Monte Carlo simulation of the full detector and offline analysis criteria. Based on 
this, we calculate 

r^cT ^T-t = (20-8 ± 0.2 ± 0.2) x lO" 3 (5) 

1 [K^ — > 7T + 7T j 

(where the first error is statistical and the second systematic) for K\ — > 7r + 7r~7 events 
with £7* > 20 MeV. The systematic uncertainty is due primarily to the effects of K s 
contamination in the vacuum beam. Assuming no contribution from interference and using 
the DE/(DE + IB) result from Table Q, we obtain the final 7r + 7r~ -normalized branching 
ratios (14.2 ± 0.2 ± 0.2) x 10" 3 for Direct Emission and (6.6 ± 0.2 ± 0.2) x 10~ 3 for Inner 
Bremsstrahlung. 

Our measured IB branching ratio is consistent with both the Q.E.D. prediction of 7.00 x 
10~ 3 and the most recent experimental result of FNAL E731 [fTOfl ; (7.3 ±0.4) x 10~ 3 , based on 
a sample of 3,136 — > 7r + 7r~7 events. The DE/(DE + IB) ratio is also in good agreement 
with the E731 result of 0.685 ±0.041, although the DE branching ratio differs from the E731 
result of (15.7±0.7) x 10~ 3 by 1.8 standard deviations. The present form-factor results differ 
significantly from that reported by E731 (ai/a 2 = —1.8 ± 0.2 GeV 2 ) but this discrepancy 
has been understood: The E731 a\ja-i form-factor was inferred |16|] from the model of Lin 



and Valencia ||, on the basis of the measured DE branching ratio, whereas our results are 
obtained directly by performing fits to the E* distribution. The underlying datasets in the 
two experiments are consistent with each other, and a reanalysis of the E731 data using our 
method yields results consistent with ours. The KTeV K° L — > 7r + 7r~7 DE form-factor result 
is also in good agreement with the result ai/a2 = —0.720 ± 0.029 GeV 2 extracted from the 
independent K\ — > 7r + 7r~e + e~ analysis from the KTeV E799 data set |4]]. 

In conclusion, we have made the first direct measurements of the — > 7r + 7r~7 Direct 
Emission form-factor, including axja-i = —0.737 ± 0.034 GeV 2 . We find no evidence for 
new CP-violating effects in the photon energy spectrum. Finally, we have made improved 
measurements of the DE ((14.2 ± 0.2 ± 0.2) x 10" 3 ) and IB ((6.6 ± 0.2 ± 0.2) x 10~ 3 ) 



K\ — > 7T + 7r _ 7 branching ratios, normalized to the K\ — >• 7r + 7r~ channel. 
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FIGURES 




FIG. 1. The decay — > ir + n~j arises primarily from the contributions of the Electric Dipole 
(El) "Inner Bremsstrahlung" (left) and Magnetic Dipole (Ml) "Direct Emission" (right) diagrams 
above. 





M,„ (GeV/c J ) nny P, 2 (GeV 2 /^) 

FIG. 2. a) distribution of candidate K® — > 7r + 7r~7 events, all other cuts applied. Arrows 

indicate final cuts at ± 8 MeV/c 2 . b) Kaon transverse momentum squared (Pj*) distribution 
of candidate — » 7r + 7r~7 events, all other cuts applied. The cut on this quantity requires events 
to be within the first bin (0.000250 GeV 2 /c 2 ) on this plot. 
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FIG. 3. K\ — ► 7r + 7r _ 7 Monte Carlo/data overlay of photon energy distribution in the center 
of mass, for the best fit result. Shown also (shaded) are the expected distributions for pure El 
Inner-Bremsstrahlung and form-factor-modified Ml Direct Emission. The "combined" Monte 
Carlo plot shown assumes these two are the only contributions to the decay. 
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FIG. 4. Ratio (points) of (IB-Subtracted) Direct-Emission data to the expectation for a 
pure Ml E* spectrum. Vertical scale is arbitrary Shown for comparison are the p-propagator 
(solid), linear (dashed), and quadratic (dotted) form-factor parametrizations. A modification to 
the pure-Mi spectrum is clearly supported by the data. 
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TABLES 



Quantity 


p— propagator 


Linear 


Quadratic 


X 2 /DOF 


38.8/27 


43.2/27 


37.6/26 


ai/a 2 


(-0.737 ± 0.026 ± 0.022) GeV 2 






r 




-1.739 ±0.062 ±0.052 


-2.93 ±0.41 ±0.34 


s 






3.31 ±1.15 ±0.96 


f 


0.683 ± 0.009 ± 0.007 


0.682 ± 0.009 ± 0.007 


0.684 ±0.011 ±0.007 



TABLE I. Summary of K\ — > it + it 7 physics results for the three form-factor parametriza- 
tions. First uncertainty is statistical, second systematic. 
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